function sk_output=sk(input_sig, a, k)
    total=length(input_sig(:,2));
    sk_output(1,2)=0;
    sk_output(2,2)=0;
    for i=2:total-1
        sk_output(i+1,2)=((k*(3-a)+2*k*k)*sk_output(i,2)-k*k*sk_output(i-1,2)+a*input_sig(i+1,2))/(1+k*(3-a)+k*k);
        sk_output(i+1,1)=input_sig(i+1,1);
    end
end